function traceSurface(vec, dx, fx, dy, fy)

[x,y]=meshgrid(dx:(fx-dx)/10:fx,dy:(fy-dy)/10:fy);
dir = [-vec(1)/vec(4) -(vec(2)/vec(4)) -(vec(3)/vec(4))];

z =  dir(1) + dir(2)*x + dir(3)*y;


%figure;
surface(x,y,z, ones(10));

end